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Abstract 

We show how experimental values of the lifetime of the Is level of pionium 
and of the difference between the energies of the 2s and 2p levels yield values 
of the elements aoc and Occ respectively of the s-wave scattering matrix for the 
two-channel (tt+tt^, vr'^vr'^) system at the tt^tt' threshold. We then develop a 
method, using energy independent hadronic potentials which reproduce the 
best available pion-pion phase shifts up to 500 MeV total energy in the cm. 
frame, for obtaining the values of the isospin invariant quantities Oq — Og 
and 2ao + Oq from Ooc and Occ respectively. We emphasise that the isospin 
invariant scattering lengths Qq and Oq universally used in the literature cannot 
be considered to be purely hadronic quantities. 
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I. INTRODUCTION 



The TT+TT" atom (pionium) has become of considerable experimental and theoretical 
interest. Such atoms have been clearly observed at Dubna and a lower limit obtained 
for the lifetime tis of the lowest (Is) level 0. A large collaboration is engaged in an 
experiment at CERN to measure with an accuracy of about 10% [0. In the proposal for 
this experiment the possibility of measuring the difference AW = W2s — W^2p between the 
energies of the 2s and 2p levels is also considered. The method was proposed by Nemenov |^ 
and its feasibility can be tested using the same experimental setup as for the measurement 

of Tis. 

A complementary approach has been suggested in Ref. p by a group at lUCF. They 
have considered the possibility of measuring the ratio of the rates for the decay of tt+tt^, 
either in ns bound states or low-lying continuum states, into tt^tt^ and 77. The known rate 
for the almost purely electromagnetic 77 decay would then yield the rate for the 7r°7r° decay. 
Ref. indicates however that a first search for a suitable reaction with tt+tt^ pairs in the 
final state has not yielded promising results. We therefore concentrate in this paper on the 
results expected to come from the CERN experiment. 

A measurement of Tis immediately yields the width of the Is level via the relation 

rf/(eV) = 0.658 212 2[ri,(fs)]-\ (1) 

The partial width for decay to 77 is very small, and can be derived from QED for bosons 
p. The result is 

Ti: = ^aVc = 7.2 X 10-^eV, (2) 

a being the fine structure constant and /ic the mass of a charged pion. Subtracting Eq. (H) 
from Eq. (||) gives the value of T^"'^", which we shall henceforth abbreviate to T. Since T is 
expected to be about 0.177 eV and will be measured with an accuracy of around 10%, no 
refined corrections to the result in Eq. (|^) are needed. 
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We expect then that values of F and, later, AW will become available from experiments 
on pionium. In Section 2 we shall go briefly through the formalism for pionium, which is 
identical with that given by Rasche and Woolcock [0 for pionic hydrogen. We shall give the 
explicit formulae for aoc(fni) and acc(fm) in terms of r(eV) and AW(eV) respectively. The 
matrix 

/ \ 

'^cc (^oc , . 

(3) 

y (^oc ^oo j 

is the s-wave scattering matrix for the two-channel (tt+tt^, tt'^tt^) system at the 7r"''7r~ thresh- 
old. The subscript c refers to the tt+tt" channel, o to the 7r°7r° channel. 

The challenge we are addressing is to connect these numerical results, derivable in princi- 
ple from experiments on pionium, with our knowledge of pion-pion scattering at low energies 
(by which we shall mean that the total energy W in the cm. frame is less than about 500 
MeV). This knowledge comes from two sources, experiments in which pion-pion pairs in the 
final states of various reactions are studied, and chiral perturbation theory (xPT), which 
gives a low-energy representation of the vrvr scattering amplitudes based on QCD. The diffi- 
culty of the problem lies in the subtle interplay of the two interactions present, the hadronic 
and the electromagnetic. The values of Oqc and a^c which come from pionium experiments 
are not due to the hadronic interaction alone: they are affected by the electromagnetic in- 
teraction. On the other hand, xPT as it has been developed so far, with the quark mass 
difference neglected and only the QCD lagrangian used, is a hadronic interaction theory 
which is charge independent (isospin invariant). Moreover, experiments above the vr+vr" 
threshold are analysed to give what are quoted as phase shifts for states of definite isospin 
((5q and 5o foi' ^=0, /=0 and 2; b\ for £=1, / =!;■ ■ ■). In spite of this, as we shall see, neither 
xPT as it is presently applied, nor the analysis of experiments, completely disentangles the 
electromagnetic interaction from the hadronic. Nevertheless it is clear that electromagnetic 
corrections need to be made, to relate the values of Ooc and obtained from pionium to 
the values of the isospin invariant s-wave scattering lengths Oq and Oq obtained from above 
threshold experiments and xPT. 



One approach has been made in a recent preprint by Maltman and Wolf , who attempt 
to include the electromagnetic interaction as an additional symmetry breaking term in xPT 
in order to calculate the corrections to the purely hadronic tttt scattering lengths. At this 
stage of their work too few details are given to fully evaluate its consequences for pionium. In 
this paper we use a potential model, considered previously by Moor, Rasche and Woolcock 
IP, to calculate the electromagnetic corrections described above. 

Moor et al. |^ attempted to calculate aoc — ^ {al — a^) by distinguishing between the 
physical situation and the purely hadronic situation in which the electromagnetic interaction 
is switched off (a = 0). In the first, the charged and neutral pions have their physical masses 
fic and fio respectively and the hadronic and electromagnetic interactions are present. In 
the second, there is only the hadronic interaction and, in accordance with current opinion 
10|-|T2[, it is assumed that all the pions have almost the same mass, which is extremely 



close to fio , and that there is practically no isospin breaking of the hadronic scattering 
amplitudes. Ref. connects these two situations by means of a potential model. For the 
purely hadronic situation, all the pions were taken to have the mass fio and only a charge 
independent hadronic potential was present. Simple square well potentials were used in a 
relativistic Schrodinger equation in order to reproduce at the energy W = 2/io the purely 
hadronic scattering lengths Oq, / = 0, 2, obtained from xPT. For the physical situation a 
two-channel s-wave relativistic Schrodinger equation was used, with physical masses and 
(relativistic) cm. momenta at a fixed energy, the potential matrix being the sum of a 
hadronic potential and a Coulomb potential for extended charge distributions. In the spirit 
of xPT it was assumed that the mass difference is the only source of symmetry breaking, so 
that the hadronic potential for the physical situation was taken to be exactly the same as the 
charge independent potential for the purely hadronic situation. The two-channel equation 
for the physical situation was then solved to obtain aoc at = 2/ic, and thus the difference 

There are two serious problems with the work in Ref. First, the purely hadronic 
scattering lengths aj] and Oq were obtained by inserting the mass fio into the xPT equations 



of Gasser and Leutwyler However some constants in those equations are fixed by using 
experimental data extrapolated to the tt+tt^ threshold W = 2fic- This is an inconsistent 
procedure. In the present state of our knowledge it is not possible to obtain purely hadronic 
scattering lengths at the purely hadronic threshold W = 2fio- One is prevented by the fact 
that, in the analysis of experiments and in the practical use of xPT, the vrvr system at the 
hadronic level is treated as an isospin invariant system with all the pions having the mass 
fic- On the basis of this assumption charge independent phase shifts are obtained from 
experiments. This information is then used, in conjunction with the xPT equations (the 
pion mass being taken as /ic), to generate s-wave phase shifts 6q and 6q for W > 2fic. Sets 
of these phases have been provided to us by Gasser . We note here that their behaviour 



at low energies is different from that in the fits of Froggatt and Petersen and Lohse et 



al. |T^. Using the expansion 

tan = + "o^c + " " " , ^ = 0, 2, (4) 

the xPT equations show clearly that > 0, < so that the curvatures have the 
same signs as the respective scattering lengths Oq. For this reason, we have reservations 



about the potentials constructed by Sander, Kuhrts and von Geramb |T6[, based on the fits 
of Refs. 

The second problem in the work of Ref. is that simple square well potentials were 
used. To obtain the values of both Qq and given by xPT, it is necessary for such potentials 
to be energy dependent. Further, Moor et al. calculated the difference between the purely 
hadronic scattering matrix at W = 2fio and the physical scattering matrix at W = 2fic 
without taking into account the variation of their potentials with energy. This particular 
problem no longer arises in the present work, where we shall calculate electromagnetic 
corrections a.t W = 2fic- Nevertheless, there remain strong reasons for desiring an energy 
independent hadronic potential matrix 



V 



had 



(5) 



If such potentials Vq can be found, the representation of the strong interaction by a potential 
matrix which remains unchanged when the electromagnetic interaction is switched on is 
placed on much firmer ground. The assumption that the tttt system at the hadronic level can 
be treated as an isospin invariant system with all the pions having the mass fic then translates 
into a two-channel Schrodinger equation for which the equations decouple in the isospin 
basis. The mass fic plays the role of a convenient reference mass and we emphasise that the 
scattering lengths Oq obtained from the potentials Vq by using a Schrodinger equation with 
this mass are isospin invariant but not purely hadronic quantities. 

There is a further question about how to construct energy independent potentials Vq . If 
the phase shifts 5q were perfectly known and the process remained elastic at all energies, it 
would be possible by inversion to construct potentials Vq . However, there are substantial 
uncertainties in the values of 6q (and electromagnetic corrections in the sense of Ref. ||17 



have not been made), there are inelastic processes at sufficiently high energy and, most 
important of all, there is no reason to believe that the hadronic interaction can be represented 
by potentials except at low energies. We therefore do not use the results of Ref. |]TB|] but 
instead find potentials of simple shape, characterised by a small number of parameters, 
which can be fixed so that the potentials reproduce the phase shifts well up to W = 500 
MeV. In fact double square wells with in all four parameters work quite satisfactorily. A 
sum of two gaussians or of two lorentzian shapes were also tried, but the fit to the phase 
shifts could not be further improved and it was more difficult to ensure that the potentials 
had no bound states. 

The double square well energy independent potentials Vq are fed into the one-channel 
Schrodinger equations 



+ 5c - /^cn'(r) j u{r) = 0, / = 0, 2, 



(6) 



and the ranges and depths of the two wells adjusted to obtain the best possible fit to the 
phase shifts 6q up to W = 500 MeV. Eq. (^) also yields the scattering lengths Oq and we 
denote by a the matrix 
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gag -t- gO-o 3 l"o 



y 3 



.ON 



l/,0 I 2 2 

3"o "T 3"o y 



(7) 



For the physical situation we have the two-channel Schrodinger equation with physical 



masses 



rf2 



l2^ + Q - 2MV^"'*(r) - 2MV^(r) 1 u(r) 



(8) 



where V*^**^ is given by Eq. (Kf) and 



Q = 


/ 

(Ic 




,2M = 


/ 


o] 


,V<"-(r) = 




"1 












f^o ) 




y » 





The relativistic cm. momenta Qc, Qo in the two channels are calculated from the same value 
of W] ^*^(r) is the Coulomb potential between extended charge distributions. The coupled 
equations (H) are integrated numerically to give the scattering matrix K as a function of 



W. The details are exactly the same as those given in Section 3 of Ref. for the coupled 
(tt^P, TT^n) system; the matrix A introduced in Eq. (65) of that paper we now denote by 
K, in accordance with the altered notation of Ref. 0. Extrapolation to = 2/ic gives the 
threshold scattering matrix a = K (2yUc) of Eq. (|^) and thus the difference a — a, which 
is the final goal of our calculation. The values of aoc and acc obtained from experiments 
on pionium can then be corrected to give values of (cq — Qq) and (20° + Oq) which can be 
compared with the values obtained from above threshold experiments and x PT. In Section 
3 we shall give further details of the calculations sketched above, and the numerical results 
for the electromagnetic corrections. 



II. DETERMINATION OF SCATTERING PARAMETERS FROM PIONIUM 

DATA 

The results in this section for the positions and widths of the pionium states of interest 
come from the formalism of Rasche and Woolcock 0], which uses analytic continuation of 
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the scattering matrix below the threshold for a channel containing two oppositely charged 
hadrons. These formal results have been confirmed by determining the position and width 
directly, by studying the behaviour of the phase shift for the (open) tc^tc^ channel. We also 
remark that the treatment of Ref. |^ is complete, and includes what are often (misleadingly) 
called hadronic effects on the wavefunction at the origin. The results of Ref. for the 
ns states of pionium, which are quite different from ours, are known to be wrong. New 



results for the Is state are given in a recent preprint by Sander and von Geramb ||I9[, but 
no comment is made on the dramatic difference between these new results and those of Ref. 
0- 

The width F for the decay of pionium in the Is state to vr^vr^ is given in Eqs. (2.10) and 
(2.13) of Ref. as 
„ 8 



B 



\l+p,{ls)eo{ls) + ■■■ \'[-Wc{l)]qXc/{^ + qXo), (9) 



where 



B = ( 2"/"c) , -Wc{l) = ^a^fic, 



pi(ls) = ^(1 +7),Reeo(ls) = ^ajB, 



a being the fine structure constant and 7 Euler's constant. The treatment of the ir'p atom 
in Ref. 0, from which the result, Eq. (^, is taken, takes account neither of the relativistic 
correction to the bound state position nor of the presence of the vacuum polarisation poten- 
tial in the calculation of F. The first is of no consequence at the level of accuracy with which 
we are concerned. A full treatment of the second has been worked out by Oades, Rasche 
and Woolcock pO| and leads to a small modification of Eq. (P). The width in the presence 
of vacuum polarisation is increased by the factor (1 + 3.08 x 10"'^); this result has been 
obtained both analytically and by a computer calculation of the width with and without the 
vacuum polarisation potential. 

Strictly speaking, the elements of the scattering matrix appearing in Eq. (|^) should be 
the elements of the K-matrix at the bound state position, not at the tt+tt^ threshold. The 
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only significant effect is on K^^] the small modification to Eq. (^) is estimated in Ref. 
and amounts to a factor (1 — 2.4 x 10~^). 

Eq. (|^) requires the values of acc and Ooo, for which we use the final numbers from Section 

3, 

a^c = 0.148 fm, a^^ = 0.053 fm. 

With the other numbers 

B = 387.489 fm, -Wc{l) = 1858.073 eV, 
Qo = 35.509 MeV, pi(ls) = 0.788608, 
he = 197.327053 MeV fm, 

we obtain the result 

aoe(fm) = 0.3802[r(eV)]i (10) 

The uncertainty in the conversion constant of Eq. (|I^) may be roughly estimated as about 
0.02%. Most of it comes from the uncertainty in the value of acc used to calculate Reeo(ls). 
Since the expected experimental accuracy in is about 5%, the uncertainty in the conver- 
sion constant is of academic interest only. 

The difference AW between the energies of the 2s and 2p levels has contributions from the 
relativistic corrections, the vacuum polarisation potential and the presence of the hadronic 
interaction. From the results of Austen and de Swart [21] the relativistic corrections to the 
positions of the 2s and 2p states are 

-37/icaVl024 = -0.014 eV(2s), 
-31/XcaV3072 = -0.004 eV(2p). 
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Using a vacuum polarisation potential modified at short distances to take account of the 
extended charge distributions of the pions, the corresponding shifts in the positions of the 
2s and 2p states are 

-0.110 eV(2s), -0.004 eV(2p). 

Thus 

AW'^^'^eV) = AW{eV) + 0.116, (11) 

where AW^"''^ is the difference due to the hadronic interaction. 

For the shifts in the positions of the 2s and 2p states due to the hadronic interaction we 
again appeal to the results of Ref . |^ . For the 2p state the shift is 

which is completely negligible. Thus AW^"''^ comes entirely from the shift in the 2s state, 
AW^'^'^^ = (-^aV) (^) (l+Pi(2s)^ + ■■■), (12) 

where 

pi(2s) = ln2 + 7 - ^ = 0.770363. 



It follows from Eq. ([TD that 

acc(fm) = -0.4168Aiy'*'^'^(eV). (13) 

From our present knowledge of a^c we expect that AW^"''^ will be close to —0.355 eV. 
The experimenters therefore face the difficult task of measuring with reasonable accuracy 
a difference AW which is about —0.47 eV, compared with the position of the n=2 pure 
Coulomb levels 464.52 eV below the tt^tt" threshold. 
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III. ELECTROMAGNETIC CORRECTIONS TO THE THRESHOLD 

SCATTERING MATRIX 



The final step is to calculate the difference a — a , where a and a are defined in Eqs. 
and (|^ respectively. As explained in Section 1, the first step is to construct energy inde- 
pendent double square well potentials, using the Schrodinger equation which reproduce 
as accurately as possible up to 500 MeV the phase shifts 6q and 6q provided by Gasser [jT^. 
There are inner and outer ranges of the potentials to be varied, and it was found that the 
best fits are obtained with ranges 



Ti = 0.3 fm, r2 = 1.5 fm. 



The full potentials which we used to obtain the final results are, to 4 significant figures. 



-6824 MeV, < r < 0.3 fm, 

+ 102.0 MeV, 0.3 f m < r < 1.5 fm, 

0, r > 1.5 fm; 



(14) 



V^{r) = +39980 MeV, < r < 0.3 fm, 

= -53.57 MeV, 0.3 f m < r < 1.5 fm, 



0, r > 1.5 fm. 



(15) 



The corresponding scattering lengths, which agree with those obtained by extrapolating 
the input phase shifts to W = 2fic, are 



a° = 0.2883 fm , al = -0.0617 fm , 



(16) 



and the matrix a of Eq. (^ is 



^ 0.1716 fm -0.1650 fm ^ 
-0.1650 fm 0.0550 fm 



(17) 
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Fig. 1 shows that the double square well potentials of Eqs. (|T4D and (0) reproduce 
the phase shifts of Ref. [|1^] remarkably well. Not surprisingly, this agreement deteriorates 
above 500 MeV but this is of no concern. The simple effective potentials clearly represent 
the hadronic interaction very well over an energy range which is quite sufficient for our 
purpose. 

The scattering lengths (Eq. (|T6|) ) come predominantly from the inner parts of the po- 
tentials in Eqs. ( |14|) and ( |T5|) and the inner radius of 0.3 fm is approximately the smallest 
radius used by Moor et al. |^ for their single square well potentials. The use of the much 
larger outer radius of 1.5 fm with the small outer potential, is a simple way of getting good 
fits to the phase shifts. The scattering lengths in Eq. ([T6|) differ from those of Ref. mainly 
because the reference mass used in the xPT equations explained in Section 1. 

The Coulomb potential V^'^(r) for extended charge distributions was taken as the poten- 
tial between uniformly charged spheres, exactly as in Ref. (last paragraph of Section 4). 
We have checked that using gaussian charge distributions instead makes no difference to the 
scattering matrix in Eq. (|TBp at the level of accuracy given. 

With the input potentials V'-'{r) just described and V*^*^*^ obtained from Eqs. (|), (|T3p 
and ([T3|), the threshold scattering matrix arising from Eq. is 



^ 0.1479 fm -0.1599 fm ^ 
-0.1599 fm 0.0533 fm j 



(18) 



The desired electromagnetic corrections are therefore 

Sec — = 0.0237 fm , doc — o-oc = —0.0051 fm, (19) 
or, with the definitions 

(Xcc (1 ~l~ ^cc)^cci ^oc (1 ~l~ ^oc)^oci 

Acc = 0.160, Aoc = 0.032. (20) 
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It is of interest to note that the contribution to the corrections from the pion mass difference 
is far greater than that from the Coulomb potentiaL The detailed numerical results are given 
in Table 1. 

On the basis of many calculations which we have made using different scattering lengths 
Oq and Oq as the starting values, we have found that the relative correction Aq^ is much less 
sensitive to these starting values than is the difference hoc — doc- However, the difference 
Occ — is less sensitive than Ace- The uncertainties in the values of Occ — dec and Aoc need 
to be estimated. Even though some fine details of the potentials have no doubt been missed, 
we are confident that the use of double square well potentials introduces a negligible error. 
We have estimated the uncertainties due to the variation of the hadronic potentials within 
limits which correspond to the uncertainties in the input phase shifts and in our fitting 
procedure. This effectively includes the uncertainties in the input scattering lengths. The 
final results are 

acc - acc = 0.024 ± 0.002 fm, A^^ = 0.032 ± 0.002, (21) 

where we have kept the number of significant digits appropriate to the size of the estimated 
errors. The development of two-loop xPT will result in slightly different phase shifts, 
but the resulting changes in dec — dec and in Aoc will lie within the uncertainities given in 
Eq. dl]). 

The values of Oqc and a^c obtained from experimental data on pionium must be corrected 
using Eq. ( pi]) , and the resulting values of hoc and dec compared with the values of Og and 
Oq obtained from a combination of vrvr scattering data at energies above W = 2/ic and the 
constraints of xPT. If this comparison should eventually reveal a significant discrepancy, it 
will then be necessary to look again at what is known about vrvr scattering above the tt+tt" 
threshold. 
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FIGURES 

FIG. 1. The 1 = and / = 2 tttt s-wave phase shifts 6q and 6q given by xPT (Ref. and 
by the double square wells of Eqs. ( p^ and (15). 
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TABLES 

TABLE I. The separate effects of the pion mass difference and of the Coulomb potential on 

the elements of the scattering matrix at the tt^tt^ threshold (all results in fm) 

(^cc ^oc ^oo 

hadronic 0.1716 -0.1650 0.0550 

with Coulomb potential only 0.1719 -0.1662 0.0561 

with mass difference only 0.1478 -0.1589 0.0522 

fully corrected 0.1479 -0.1599 0.0533 
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